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ABSTRACT 


Some  of  the  theoretical  and  experimental  implications  of  isotropic 
incremental  elastic-plastic  constitutive  models  formulated  with  constant 
values  for  the  shear  modulus  G  are  examined.  Two  types  of  constant  G 
models  are  considered,  i.e.  a  single  constant  value  of  G  for  loading  and 
unloading  and  two  separate  values  of  G  ,  one  for  loading  and  one  for  un¬ 
loading.  These  material  property  specifications  are  successively  coupled 
with  four  increasingly  more  realistic  idealized  stress-strain  relations 
representing  constrained  modulus  M  functions  obtained  for  a  state  of 
uniaxial  strain  and  a  von  Mises-type  limiting  shear  envelope  characteristic 
of  that  specified  fc  the  classical  Prandtl-Reuss  material.  The  resulting 
models  a^e  used  to  calculate  mean  normal  stress  or  pressure  versus  volu¬ 
metric  strain  and  uniaxial  strain  principal  stress  difference  versus  pres¬ 
sure  relations  for  qualitative  comparison  with  observed  test  phenomena. 

The  results  from  this  simple  examination  of  several  highly  idealized 
material  descriptions  indicate  that  many  of  the  behavior  characteristics 
of  real  earth  materials  often  recorded  during  laboratory  tests  can  be  simu¬ 
lated,  at  least  qualitatively,  by  elastic-plastic  constitutive  models  of 
the  constant  G  type  by  observing  certain  restrictions  on  the  material 
property  parameters  G  and  M  .  Other  characteristic  phenomena,  however, 
cannot  be  mirrored  with  these  types  of  models  under  any  circumstances. 

The  constitutive  equations  of  a  classical  Prandtl-Reuss  material  and 
the  behavior  of  this  ideal  model  for  conditions  of  uniaxial  strain  are  in¬ 
cluded  in  an  appendix  for  reference  and  informational  purposes. 
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PREFACE 


The  work  described  herein  is  part  of  a  continuing  research  effort 
being  conducted  by  personnel  of  the  Soils  Division,  U,.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  aimed  at  developing  physically  re¬ 
alistic  and  theoretically  satisfying  constitutive  models  for  a  broad 
range  of  soil  and  rock  materials  for  use  as  input  to  free- field  ground 
shock  calculation  computer  codes.  This  particular  study,  which  deals  in 
general  terms  with  relatively  simple  idealizations  of  material  behavior, 
developed  as  an  adjunct  to  an  investigation  begun  by  CPT  J.  A.  Spitz- 
nagel,  CE,  into  the  use  of  constant  shear  modulus- type  constitutive 
models  to  represent  the  actual  behavior  observed  in  the  WES  Soil  Dynamics 
Laboratory  during  numerous  tests  on  soils  and  weak  sedimentary  rocks . 

The  theoretical  analyses  of  these  idealized  models  were  performed  during 
the  period  February- April  1971  by  Dr.  Behzad  Rohani;  this  paper  was 
written  by  Dr.  Rohani  and  Dr.  J.  G.  Jackson,  Jr. 

The  project  is  sponsored  by  the  Defense  Nuclear  Agency  as  part 
of  their  Nuclear  Weapons  Effects  Research  Subtask  SB209,  "Propagation 
of  Ground  Shock  Through  Earth  Media."  The  work  is  being  performed 
under  the  genb^ai  direction  of  Messrs.  J.  P.  Sale,  R.  G.  Ahlvin,  and 
R.  W.  Cunny,  of  the  Soils  Division,  WES.  The  Director  of  WES  during  the 
preparation  and  publication  of  this  report  was  COL  Ernest  D.  Peixotto, 

CE.  Technical  Director  was  Mr.  Fred  R.  Brown. 
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CHAPTER  1 


INTRODUCTION 

1.1  BACKGROUND 

The  constitutive  equations  currently  used  to  represent  various  types 
of  earth  media  in  two-dimensional  ground  shock  calculation  codes  are 
based  on  the  assumption  of  isotropic  incremental  elastic- plastic  material 
behavior.  In  the  elastic  range,  material  behavior  is  described  by  the 
following  incremental  constitutive  relation  that  effectively  separates 
the  deformation  into  hydrostatic  and  deviatoric  components. ^ 


ct.  .  =  Ke, ,  5.  .  t  2G(e.  .  -  k  e,,8.  ..) 

10  kk  io  v  10  3  kk  io 


(1.1) 


where  a. .  =  components  of  the  total  stress  increment  tensor 
id 

=  ■—  a, .  6 .  .  +  s.  . 

3  kk  io  io  . 

K  =  bulk  modulus  =  -r^~ 
ekk 

=  sum  of  the  incremental  normal  strains  =  incremental  change 
in  volumetric  strain 


5 .  .  =  Kronecker  delta 

s 

G  =  shear  modulus  =  t, 

-  components  of  the  total  strain  increment  tensor 

1  *  * 

=  -  t,6.  ,  +  e.  . 

3  kk  io  ij 

t  • 

=  sum  of  the  incremental  normal  stresses  =  3P 

• 

sio  =  comPonen^s  incremental  stress  deviator  tensor 

e^  =  components  of  the  incremental  strain  deviator  tensor 


Elastic  range  as  used  in  this  paper  does  not  imply  linear  reversible 
behavior  of  the  material,  but  rather  refers  to  all  stress- strain  re¬ 
sponses  for  deviatoric  stress  states  less  than  that  defined  by  an  ex¬ 
plicit  plastic  yield  function. 
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llonlinear  behavior  is  incorporated  by  making  K  and  G  functions  of 
stress  and/or  strain  invariants.  Hysteretic  behavior  is  taken  into  ac¬ 
count  usin'*  different  values  of  K  and  G  on  loading  and  unloading. 

A  pseudoplastic  behavior  is  tlms  built  into  the  incremental  constitutive 
relation  (Equation  l.l)  in  the  form  of  permanent  compaction  prior  to 
yield . 

Classical  plasticity  is  incorporated  into  the  model  by  specifying  a 
yield  condition  that  effectively  serves  to  limit  the  maximum  shear 
strength  of  the  material.  The  yield  surface  is  usually  a  combination  of 
the  li-ager- Trucker  and  von  Mises  types,  often  with  a  functional  form  that 
ensures  a  smooth  transition  from  the  former  to  the  latter  with  increasing 
pressure  (References  1  and  2) . 

With  zero  work  hardening,  the  yield  function  is  specified  by  a  rela¬ 
tion  of  the  form 


f<>i\  =  f(p) 

max 


C  1.2) 


where  JA  =  second  invariant  of  the  stress  deviator  tensor  -  s.  .s.  . 

2  2  xj  ij 

p  =  mean  normal  stress  or  pressure  = 

sij  =  comPonen't's  °-  the  stress  deviator  tensor 

From  within  this  generalized  description  of  elastic-plastic  material 
behavior,  many  specific  descriptions  are  available,  depending  upon  the 
various  restrictions  placed  on  the  constitutive  property  parameters  K  , 

One  of  the  more  common  versions  involves  consider- 


0  > 


and 


wl)max 

able  flexibility  in  describing  volume  change  characteristics  by  permit- 
ing  K  to  be  defined  by  nonlinear  polynomial  functions  while  restricting 
the  description  of  shear  behavior  to  constant  values  for  G  and 
The  significance  of  such  restrictions  can  perhaps  best  be  illustrated  by 
examining  the  behavior  of  these  constant  shear  modulus  models  when  sub¬ 
jected  to  the  states  of  stress  most  often  employed  in  the  laboratory  to 
experimentally  determine  the  stress- strain  response  of  actual  earth 
materials . 


The  testing  technique  most  frequently  used  for  studying  stress- strain 
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response  in  the  laboratory  in  conjunction  with  ground  shock  problems  in¬ 
volves  confiring  undisturbed  specimens  so  that  they  deform  in  an  undrained 
state  of  uniaxial  strain  under  an  applied  transient  stress  (Reference  3) . 
The  test  results  are  usually  illustrated  through  a  plot  of  vertical  stress 

versus  vertical  strain;  the  slopes  of  such  plots  define  the  constrained 
2 

modulus  M  ,  i.e. 


M  = 


\,e  s/uniaxial 
strain 


(1-3) 


Triaxial  shear  tests  are  normally  used  to  determine  relationships  between 
JfiT)  or  maximum  principal  stress  difference  (a  -  a )  ,  and  p  . 
For  a  given  material,  a  unique  path  of  (a  -  a  )  versus  p  is  also  de- 
scribed  by  a  material  undergoing  uniaxial  strain  loading  and  unloading; 
ohe  slopes  along  a  stress  path  of  this  type  define  the  ratio  of  shear 
modulus  to  bulk  modulus ,  i.e. 


'■uniaxial 

strain 


2G 

K 


(1.4) 


Thus,  as  an  aid  in  -understanding  the  practical  significance  of  constant 
G  constitutive  models,  it  is  helpful  to  examine,  for  a  given  uniaxial 
strain  o_  versus  c  relation  and  a  given  triaxial  shear  envelope,  the 
characteristic  features  of  the  corresponding  p  versus  €  and  uniaxial 
strain  (a  -  aj)  versus  p  relationships  imposed  by  the  further  specifi- 
cation  of  a  constant  value  for  G  . 


1.2  PURPOSE  AMD  SCOPE 

The  purpose  of  this  paper  is  to  examine  some  of  the  theoretical  and 


The  laboratory  tests  currently  used  for  soil  property  investigations 
are  axially  symmetric  so  that  the  specimens  can  be  analyzed  in  terms 
of  coordinates  s  and  r  denoting  vertical  and  radial  directions, 
respectively • 
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experimental  implications  of  isotropic  incremental  elastic-plastic  consti¬ 
tutive  models  formulated  with  constant  values  for  the  shear  modulus  G  . 

Two  types  of  constant  shear  modulus  models  will  be  considered.  The  first 
model,  in  which  one  single  constant  value  of  G  is  used  for  loading  and 
unloading,  is  referred  to  as  a  "linear  elastic  G  model."  The  second, 
in  which  two  separate  constant  values  of  G  are  specified,  one  for  load¬ 
ing  and  one  for  unloading,  is  called  the  "linear  hysteretie  G  model." 

Shear  stress- shear  strain  diagrams  for  both  models  are  given  in  Figure  1.1. 

Four  idealized  uniaxial  strain  a  versus  e  relations  are  used 

z  z 

in  conjunction  with  the  two  constant  G  relations.  The  first  represents 
the  simplest  idealisation,  i.e.  that  defined  by  a  linear  elastic  con¬ 
strained  modulus .  The  second  is  in  the  linear  hysteretie  form  defined  by 
separate  constant  values  for  M  ,  i.e.  one  for  loading  and  one  for  unload¬ 
ing.  The  third  represents  a  bilinear  hysteretie  constrained  modulus  de¬ 
scription  and  the  fourth  a  general  nonlinear  hysteretie  description.  The 
four  relations  are  diagramed  in  Figure  1.2. 

It  is  assumed  that  plastic  behavior  is  governed  by  a  von  Mises-type 
yield  criterion  as  postulated  for  the  classical  Prandtl-Reuss  material. 
Chapter  2  presents  an  analysis  of  six  models,  developed  by  successively 
coupling  the  two  constant  G  descriptions  (Figure  l.l)  with  the  M  de¬ 
scriptions  outlined  in  Figures  1.2a,  b,  and  d,  subjected  to  maximum  devi- 
atoric  stress  levels  insufficient  to  invoke  the  yield  criterion..  Two 
cases  in  which  the  deviatoric  stress  paths  reach  the  plastic  yield  sur¬ 
face  are  examined  in  Chapter  35  these  cases  are  defined  by  successively 
coupling  the  linear  elastic  constant  G  descriptions  with  the  two  M 
descriptions  given  in  Figures  1.2c  and  d. 

For  reference  and  informational  purposes,  the  mathematical  constitu¬ 
tive  equations  for  a  linear  elastic-plastic  model  of  an  ideal  Prandtl- 
Reuss  material  are  presented  and  its  uniaxial  strain  behavior  is  examined 
in  Appendix  A. 
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c.  BILINEAR  HYSTERETIC  M 


J.  NONLINEAR  HYSTERETIC 


Figu->.  1.2  Vertical  stress  versus  vertical  strain  curves  for 
conditions  of  uniaxial  strain. 


CHAPTER  2 


? 


CONSTANT  SHEAR  MODULUS  CONSTITUTIVE  MODELS 
FOR  STRESS  STATES  BELOW  YIELD 


The  basic  constitutive  relationship  in  the  elastic  range  is  Equa¬ 
tion  1.1,  which  contains  two  material  constants  (or  functions)  K  and 
G  .  Equation  !>  .1  can  also  be  written  in  terms  of  constrained  modulus  M 
and  shear  modulus  G  by  utilizing  Equation  1.3,  i.e. 


«lj  ‘  MVij  + 


(2.1) 


The  procedure  adopted  in  this  chapter  to  study  the  behavior  of  the  two 
constant  G  models  for  subyield  stress  states  (elastic  range)  is  to  as¬ 
sume  that  M  and  G  are  known  for  each  model  and  then  substitute  in 
Equation  2.1  to  calculate  the  pres sure- volumetric  strain  or  K  relation 
and  the  uniaxial  strain  path  of  principal  stress  difference  versus  pres¬ 
sure.  The  calculated  behavior  is  then  discussed  in  terms  of  the  general 
behavior  observed  during  tests  on  actual  earth  materials. 

If  and  G^  denote  the  constrained  and  shear  moduli,  respec¬ 
tively,  for  virgin  loading  and  My  and  Gu  denote  the  corresponding 
values  for  unloading,  the  expressions  for  the  pressure-volumetric  strain 
relation  and  the  stress  path  associated  with  the  uniaxial  strain  test 
become 


Virgin 

loading 


*  '  (Mi  - 1  °e)^  -  V, 


(az  -  ar)  » 


2G 


2G 


kM,  -  i  G. 
i  3  h 


P  -  -v r*  P 


Unloading 
and 
reloading 


P  =  /m  -  ^  G  )  e  =  K  e 
\  u  3  u/  z  us 


(a  -  o )  = 

c*  1 


2G. 


u 

T 


2G. 


«u  -  1  % 


u  • 

P  =  —  P 
*u 


(2.2a) 


(2.2b) 


(2.2c) 


(2. 2d) 
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The  values  of  loading  and  ■unloading  shear  moduli  are  constant  throughout 

the  calculations,  whereas  the  values  of  loading  and  unloading  constrained 

moduli  could  vary  with  stress  level,  in  which  case  the  construction  of 

the  model  must  be  carried  out  incrementally  by  solving  for  compatible 

values  of  (a  -  o„)  and  p  corresponding  to  each  pair  of  c  .e,,  values. 
Z  x  z  z 

The  pertinent  incremental  relations  are 

<az  ~  °r ^  =  2G/z  for  l°adinS  (2.3a) 


-  a  )  ?  2G  e 
r'  u  z 


for  unloading- reloading 


(2.3b) 


P  =  °z  ‘  |  (oz  "  crr)  (2.3c) 

The  total  quantities  needed  for  plotting  pertinent  stress- strain  and 
stress  path  relations  are  then  obtained  from  the  following  expressions: 


Ni =  Nt.x  +  **  ('dM) 

(%  -  C*)^*'*  (2-4b) 


(2.4c) 


(p)i  -  (p)^!  +  P  (2.4d) 

where  i  =  1  ,  m  .  The  calculation  is  performed  in  each  increment 
of  loading  or  unloading  as  if  the  material  were  elastic. 

Six  different  models  are  considered  in  this  chapter  using  various 
combinations  of  constrained  and  shear  moduli  as  outlined  in  Section  1.2. 
The  models  are  presented  below  in  the  order  of  increasing  complexity. 
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2.1  LINEAR  ELASTIC  G  MODEL  FOR  MATERIALS  WITH  A  LINEAR  ELASTIC  M 


This  model  corresponds  to  a  linear  elastic  material  and  is  obviously 
the  simplest  possible  constitutive  model  available.  It  is  presented  here 
only  to  demonstrate  the  methodology  used  in  this  chapter.  Since  both  the 
shear  modulus  and  the  constrained  modulus  are  elastic,  it  implies  that 
M.  =  M  =  M  and  G  =  G  =  G  so  that  the  pres sure- volumetric  strain  re- 
lation  and  the  stress  path  associated  with  a  state  of  uniaxial  strain  for 
this  model  become 


■  (m  - 1  °K 


(2.5a) 


(os  " 


CTr}  = 


(2.5b) 


From  thermodynamic  considerations,  it  is  necessary  to  impose  the  follow¬ 
ing  inequalities  on  the  shear  modulus  G  and  on  Poisson’s  ratio  v 


G  >  0 


(2.6a) 

(2.6b) 


in  order  to  insure  the  positive  definite  character  of  the  quadratic  form 
of  the  elastic  strain  energy  function  (Reference  4).  Negative  values  of 
v  ,  however,  have  not  been  found  experimentally  for  any  isotropic  elastic 
material.  Thus,  from  purely  phenomenological  considerations,  the  re¬ 
striction  on  v  should  be  revised  to  limit  it  to  positive  values,  i.e. 


0  <  v  <2 


(2.7) 


The  relationship  between  v  ,  M  ,  and  G  is  given  from  elastic  theory 
as 


=  M  ~  2G 
V  2(M  -  G) 


(2.8) 
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The  pertinent  stress- s>  /rain  and  stress  path  plots  for  this  simple  elastic 
model  are  shown  schematically  in  Figure  2.1. 


2.2  LINEAR  HYSTERETIC  G  MODEL  FOR  MATERIALS  WITH  A  LINEAR  ELASTIC  M 

Figures  2.2a  and  b  depict  the  vertical  stress- strain  and  the  prin¬ 
cipal  stress  difference- strain  difference  relations,  respectively,  for 
this  model  under  uniaxial  strain  conditions.  Since  the  constrained  modu¬ 
lus  is  linear  elastic  (M  =  M  =  M)  and  the  shear  modulus  is  linear  hys- 

A/ 

teretic  (G  >  G.),  the  pressure- volumetric  strain  and  the  principal  stress 
difference-pressure  relations  for  the  model  are 


(2.9a) 

(2.9b) 


Equations  2.9a  through  d  are  plotted  schematically  in  Figures  2.2c  and  d. 
The  material  is  first  loaded  (in  uniaxial  strain  configuration)  to  a 
stress  level  corresponding  to  Point  1  and  then  unloaded  to  a  stress  level 
corresponding  to  Point  2. 

Figure  2.2c  demonstrates  an  important  aspect  of  the  constant  shear 
modulus  model  when  there  is  no  hysteresis  in  the  i  niaxial  strain  curve 
(M  =  M  ).  Since  G  >  G.  ,  the  unloading  slope  of  the  pressure- 

JL  U  U  a> 

volumetric  strain  curve  is  smaller  than  the  loading  slope,  i.e. 

(M  '  1  °u)  <  (“  •  3  ai)  (2-10) 


1 6 


which  results  in  an  energy- generating  loop  in  the  pressure-volumetric 
strain  curve.  This  is  in  direct  opposition  to  the  available  experimental 
evidence,  which  indicates  that  a  true  hydrostat  would  not  show  an  energy¬ 
generating  loop  and  the  model  should  not  produce  one  in  the  constructed 
pressure- volumetric  strain  curve. 

It  should  be  pointed  out  that  the  values  of  G .  and  G  are  not 

AS  ** 

arbitrary.  In  uniaxial  strain  tests,  it  has  generally  been  observed  that 
incremental  increases  in  a  are  accompanied  by  increases  in  a  and 

Z  A 

e  ,  that  incremental  decreases  in  a  result  in  decreases  in  a  and 

Z  Z  a 

e  ,  and  that  neither  a  nor  e  exhibits  tensile  values  while  sub- 
z  ’  r  z 

jected  to  compressive  values  of  a  .  Thus,  M  and  K  remain  positive 

z 

at  all  times.  This  condition  restricts  both  G^  and  G^  by  requiring 
that 


3G£  <M 


1Gu<M 


(2.11a) 

(2.11b) 


Equations  2.11a  and  b,  however,  must  be  replaced  by 


2G^  <  M  (2.11c) 

2CU  <  M  (2. lid) 

in  order  to  exclude  negative  values  of  Poisson's  ratio  (Equations  2.7  and 
2.8).  Therefore,  the  maximum  values'  of  G.  and  G  are  fixed  by  the 
value  of  constrained  modulus  M  according  to  Equations  2.11c  and  d. 

2.3  LINEAR  ELASTIC  G  MODEL  FOR  MATERIALS  WITH  A  LINEAR  HYSTERETIC  M 

Figures  2.3a  and  b  depict  the  vertical  stress- strain  and  the  prin¬ 
cipal  stress  difference- strain  difference  relations,  respectively,  for 
this  model  under  uniaxial  strain  conditions.  Since  the  shear  modulus  is 
linear  elastic  (G^  =  G^  -  G)  and  the  constrained  modulus  is  linear  hyster- 
etic  (Mu  >  M^) ,  the  pressure-volumetric  strain  and  the  principal  stress 
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difference-pressure  relations  for  the  model  are 


\ irgin 
loading 


(2.12a) 


-  o  )  =  J - SSjj — N  p  (2.12b) 

V*1  -  3  V 

(Mu'3g)%  <2-12=> 

-  %'> - ( — V  <2-i2a) 

W  -  s  v 

Equations  2.12a  through  d  are  plotted  schema, tically  in  Figures  2.3c  and  d. 

Unlike  the  previous  model,  there  is  no  energy- generating  loop  in  the 
pressure- volumetric  strain  curve  associated  with  this  model  (Figure  2.3c). 
This  is  due  to  the  fact  that  M  >  M.  ,  which  results  in  an  unloading 

U  X/ 

slope  for  the  pressure-volumetric  strain  curve  that  is  larger  than  the 
loading  slope,  i.e. 

Thus,  utilizing  a  hysteretic  constrained  modulus  in  a  constant  G  consti¬ 
tutive  model  results  in  a  pres sure- volumetric  strain  curve  that  is  quali¬ 
tatively  in  agreement  with  observed  behavior.  However,  the  constructed 
■stress  path  (Figure  2.3d)  is  contrary  to  most  of  the  existing  experimental 
data,  which  generally  indicate  that  the  unloading  stress  path  should  fall 
below  the  loading  path.  In  view  of  Equations  2.12b,  2.12d,  and  2.13,  it 
is  clear  that  unloading  stress  paths  constructed  with  this  type  model  will 
always  fall  above  the  loading  paths . 

As  discussed  for  the  previous  model,  the  value  of  G  for  this  model 
is  also  restricted  by  the  inequality 

2G  <  (2.14) 


Unloading 
and 

reloading  |  ^ 
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2.4  LINEAR  HYSTERETIC  G  MODEL  FOR  MATERIALS  WITH  A  LINEAR 

HYSTERETIC  M 

Figures  2.4a  and  b  depict  the  vertical  stress- strain  and  the  prin¬ 
cipal  stress  difference- strain  difference  relations,  respectively,  for 
this  model  under  uniaxial  strain  conditions.  Both  the  constrained  and 
shear  moduli  are  linear  hysteretic,  i.e.  >  M^  and  Gy  >  G^  .  The 
pressure-volumetric  strain  relation  and  the  principal  stress  difference- 
pressure  relation  for  the  model  are  shown  schematically  in  Figures  2.4c 
and  d.  Tne  material  is  first  loaded  to  a  stress  level  corresponding  to 
Point  1  and  then  unloaded  to  a  stress  level  corresponding  to  Point  2 
(Figures  2.4a  and  b).  The  constructed  pressure- volumetric  strain  (Fig¬ 
ure  2.4c)  and  principal  stress  difference-pressure  relations  (Figure  2.4d), 
however,  may  follow  Paths  1-A  or  1-B  upon  unloading  depending  on  the  rela¬ 
tive  values  of  M^  ,  Mu  ,  G£  ,  and  Gu  .  If  (mu  -  ^  Gu)  <  (m^  -  ^  G^)  , 
the  pressure-volumetric  strain  curve  will  unload  along  the  Path  1-A  re¬ 
sulting  in  an  energy- generating  loop,  which  is  contrary  to  existing  ex¬ 
perimental  data,  as  was  pointed  out  in  Section  2.2.  If,  on  the  other 
hand,  ^M^  -  G^  >  ^M^  -  ^  G^  ,  the  pressure-volumetric  strain  curve 
will  unload  along  the  Path  1-B,  which  is  in  general  agreement  with  the 
available  experimental  data.  The  stress  difference-pressure  relation 
will  tinload  along  Path  1-A,  which  is  contrary  to  most  of  the  experimental 

data,  if  G_J G^  <  G^/^M^  -  G^  and  will  unload  along  Path  1-B 

if  ^  G.a^  >  g/m£  -  ~  G^  .  Thus ,  in  order  for  both  the 

pres sure- vblumetric  strain  relation  and  the  principal  stress  difference- 
pressure  relation  to  unload  along  Path  1-B,  the  following  inequalities 
must  be  satisfied  for  the  model 


(\-V  >3(gu 

v  3s 


(2.15a) 


(2.15b) 


The  values  of  G^  and  Gy  for  this  model  are  also  restricted  as 
previously  noted  b;y  the  inequalities 


19 


“A  «  M* 

(2.16a) 

“u  <  t 

(2.16b) 

2.5  LINEAR  ELASTIC  G  MODEL  FOR  MATERIALS  WITH  A  NONLINEAR 

KYSTERETIC  M 

Figures  2.5a  and  b  depict  the  vertical  stress- strain  and  the  prin¬ 
cipal  stress  difference- strain  difference  relations,  respectively,  for 
this  model  under  uniaxial  strain  conditions.  The  pressure- volumetric 
strain  and  the  principal  stress  difference-pressure  relations  for  the 
model  are  plotted  schematically  in  Figures  2.5c  and  d.  This  model  is 
the  nonlinear  version  of  the  model  presented  in  Section  2.3  and  Figure 
2.3*  Although  the  principal  features  of  the  models  are  the  same,  there 
is  an  additional  feature  associated  with  the  nonlinear  model  that  is  not 
evident  in  a  linear  model.  Since  the  slope  of  the  stress  path  varies  in¬ 
versely  with  the  constrained  modulus,  it  follows  that  the  stress  path  has 
a  curvature  that  is  opposite  to  that  in  the  corresponding  portion  of  the 
uniaxial  strain  curve.  This  behavior  is  illustrated  in  Figure  2.5  by  the 
dashed  line  A-B.  The  experimental  data  obtained  to  date,  however,  indi¬ 
cate  that  the  stress  path  is  generally  continuously  concave  to  the  pres¬ 
sure  axis  throughout  the  loading  cycle  (such  as  Path  A-l  in  Figure  2.5d) 
regardless  of  the  shape  of  the  vertical  stress-strain  curve. 

2.6  LINEAR  HYSTEEETIC  G  MODEL  FOR  MATERIALS  WITH  j  NONLINEAR 

IIYSTERETIC  M 

Figures  2.6a  and  b  depict  the  vertical  stress-strain  and  the  prin¬ 
cipal  stress  difference- strain  difference  relations,  respectively.  This 
model  is  the  nonlinear  version  of  the  model  presented  in  Section  2.4  and 
Figure  2.4.  The  constructed  pressure-volumetric  strain  and  the  principal 
stress  difference-pressure  relations  for  the  model  are  plotted  schemati¬ 
cally  in  Figures  2.6c  and  d.  As  discussed  in  Section  2.4,  the  constructed 
pressure- volumetric  strain  relation  (Figure  2.6c)  and  the  principal  stress 
difference-pressure  relation  (Figure  2.6d)  may  follow  either  Path  1-A  or 
Path  1-B  upon  unloading ,  depending  on  the  relative  values  of  M  ,  M  , 

G^  ,  and  C»u  . 
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In  order  for  both  the  pressure-volumetric  strain  relation  and  the 
principal  stress  difference-pressure  relation  to  unload  along  Path  1-B, 
the  inequalities  given  by  Equation  2.15  must  be  satisfied  for  this  model 
at  any  stress  level.  To  prevent  the  occurrence  of  negative  v  values, 
the  values  of  G^  arid  Gu  must  also  be  restricted  by  the  inequalities 
given  by  Equation  2.16  for  any  stress  level. 
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a.  VERTICAL  STRESS-STRAIN 
RELATION 


b.  PRINCIPAL  STRESS  DIFFERENCE- 
STRAIN  DIFFERENCE  RE L ATION 


<°z 


c.  PRESSURE-VOL'JMETRIC 
STRAIN  RELATION 


i.  PRINCIPAL  STRESS^DIFFERENCE- 
PRESSURE  RELATION  (STRESS 
PATH) 


Figure  2.1  Behavior  of  linear  elastic  material  under 
conditions  of  uniaxial  strain. 
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VERTICAL  STRESS-STRAIN 
RELATION 


b.  PRINCIPAL  STRESS  DIFFERENCE' 
STRAIN  DIFFERENCE  RELATION 


7 

z 


c.  PRESSURE-VOLUMETRIC  i.  PRINCIPAL  STRESS  DIFFERENCE 

STRAIN  RELATION  PRESSURE  RELATION  (STRESS 

PATH) 

Figure  2.3  Behavior  of  linear  elastic  shear  modulus  model  for  materia] 
with  a  linear  hysteretic  constrained  modulus  under  conditions  of  uni¬ 
axial  strain. 


■  VERTICAL  STRESS-STRAIN 
RELATION 


b.  PRINCIPAL  STRESS  DIFFERENCE' 
STRAIN  DIFFERENCE  RELATION 


Figure  2.h  Behavior  of  linear  hysteretio  shear  modulus  model  for 
materials  with  a  linear  hysteretic  constrained  modulus  under  con¬ 
ditions  of  uniaxial  strain. 
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o.  VERTICAL  STRESS-STRAIN 
RELATION 
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b.  PRINCIPAL  STRESS  DIFFERENCE 
STRAIN  DIFFERENCE  RELATION 
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c.  PRESSURE-VOLUMETRIC 
STRAIN  RELATION 


d.  PRINCIPAL  STRESS  DIFFERENCE 
PRESSURE  RELATION  (STRESS 
PATH) 


Figure  2.5  Behavior  of  linear  elastic  shear  modulus  model  for 
materials  -with  a  nonlinear  hysteretic  constrained  modulus  under 
conditions  of  uniaxial  strain. 
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CHAPTER  3 


CONSTANT  SHEAR  MODULUS  CONSTITUTIVE  MODELS  FOR 
STRESS  STATES  REACHING  YIELD 


The  procedure  used  to  construct  the  various  constant  G  models  dis¬ 
cussed  in  this  chapter  is  basically  the  same  as  that  described  in  Chap¬ 
ter  2.  The  distinction  between  the  two  procedures  becomes  apparent  only 
when  a  stress  state  satisfying  the  von  Mises  yield  condition  =  k)  is 
obtained.  When  such  a  stress  state  is  reached,  the  value  of  shear  modu¬ 
lus  becomes  effectively  equal  to  zero  and  the  constrained  and  bulk  moduli 
become  equal.  Two  different  models  are  considered  in  this  chapter  as 
outlined  in  Section  1.2. 


3.1  LINEAR  ELASTIC  G  MODEL  FOR  MATERIALS  WITH  A  BILINEAR 
HYSTERETIC  M 

Figure  3.1a  depicts  the  vertical  stress-strain  relation  for  this 
model  under  uniaxial  strain  conditions.  The  loading  portion  of  the 
stress-strain  curve  is  approximated  with  two  straight  lines  having  slopes 

t.221 

M7  and  M.  where  M.  >  M.  .  The  unloading  part  of  the  curve  is  also 

Xt  A/  Xj  Xj  ^  g 

represented  by  two  straight  lines  with  slopes  m  and  M^  where 


<  >  Mu 


v  "  *  u  u 

The  principal  stress  difference- strain  difference  relation  for 


the  model  is  shown  in  Figure  3.1b.  The  material  first  behaves  as  a  lin¬ 
ear  elastic  G  model  with  a  linear  hysteretic  constrained  modulus  (Sec¬ 
tion  2.3)  until  a  stress  state  corresponding  to  the  von  Mises  yield  limit 
k  is  reached.  This  corresponds  to  Point  1  in  Figure  3*1*  Physically, 
Point  1  is  often  associated  with  the  forcing  of  entrapped  air  into  solu¬ 
tion  and  the  subsequent  saturation  of  the  soil  specimen,  which  results  in 
a  stiffer  oehavior  in  uniaxial  strain  or  hydrostatic  compression  and  a 
weaker  behavior  in  shear. 

Upon  yielding,  the  material  behaves  as  though  it  were  a  constrained 

fluid  exhibiting  no  shear  resistance.  The  slope  of  the  vertical  stress- 

2 

strain  curve  breaks  at  Point  1  and  becomes  equal  to  .  The  slope  of 
the  principal  stress  difference- strain  difference  curve  also  breaks  at 
Point  1  and  becomes  equal  to  zero .  The  stress  difference  and  vertical 
stress  magnitudes  at  Point  1  are,  respectively,  V§  k  and  V3  kM^/2G  . 
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During  unloading  from  a  stress  state  corresponding  to  Point  2  (Figures  3*la 

and  b),  the  material  behaves  as  a  linear  elastic  solid  in  shear  (G.  =  G 

TP  ** 

a  G),  but  since  M  >  W.  ,  it  exhibits  hysteretic  effects  in  the  hydi^ostatic 

“  *  1  4  2 

part  of  its  deformation  unless  -  —  G  =  .  When  the  lower  yield  sur¬ 

face  corresponding  to  Point  3  is  reached,  the  material  again  flows  plasti¬ 
cally.  The  slope  of  the  vertical  stress- strain  curve  breaks  at  Point  3  and 
2 

becomes  equal  to  M  and  the  slope  of  the  principal  stress  difference- 

U  2  1 
strain  difference  curve  becomes  equal  to  zero .  Since  M  >  M .  ,  hysteretic 

U  AJ 

volume  change  continues  to  increase  with  subsequent  further  unloading. 

Equation  2.1  can  now  be  utilized  to  calculate  +he  pressure- volumetric 
strain  and  the  stress  difference-pressure  relations  associated  with  a  state 
of  uniaxial  strain.  During  loading  from  Point  0  to  Point  1  (Figure  3.1), 
the  pressure-volumetric  strain  and  the  stress  difference-pressure  relations 
become 


Equations  3.1a  and  b  apply  only  up  to  Point  1,  which  corresponds  to 


(3-la) 

(3-lb) 


(az  -  ar)  -  V3k 


(3-2a) 

(3-2b) 


During  loading  from  Point  1  to  a  stress  state  corresponding  to  Point  2,  the 
material  flows  plastically  and  the  shear-  modulus  becomes  effectively  equal 
to  zero.  Accordingly,  the  pressure- volumetric  strain  relation  becomes 

(3-3) 

During  unloading  from  Point  2  to  a  stress  state  corresponding  to  Point  3, 
the  material  behaves  elastically  in  shear  resulting  in  the  following  stress 
difference-pressure  relation 


(3-4) 
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At  Point  3)  the  yield  surface  is  intersected  again,  this  time  on  its 
"lower"  side  so  that 


(°z  '  °r^  =  k 


(3-5) 


The  pres sure- volumetric  strain  curve,  however,  may  unload  along  Paths  2-A 
or  2-B  as  shown  in  Figure  3.1c,  depending  on  the  relative  values  of  the 
constrained  and  shear  moduli  used  for  the  model.  The  mathematical  expres¬ 
sion  describing  the  pressure- volumetric  strain  relation  for  the  stress 
range  between  Points  2  and  3  is  given  by 

i  -  (<  -  3  G)h  (3'6) 

If  (m'L  -  4  g)  <  ,  the  uressure- volumetric  strain  curve  will  unload 

\  u  3  /  1 

along  Path  2-A,  resulting  in  an  energy- generating  loop,  which  contradicts 

/  1  4  \  2 

available  experimental  data.  If,  on  the  other  hand,  *  —  GJ  > 

the  pressure-volumetric  strain  curve  will  unload  along  Path  2-B,  which  is 
in  general  agreement  with  current  experimental  data.  If  unloading  is  con¬ 
tinued  beyond  Point  3,  the  material  flows  plastically  again  and  the 
pres sure- volumetric  strain  relation  becomes 

i  =  (3.7) 


Equations  3*1  through  3.7  completely  describe  the  behavior  of  the 
model  under  conditions  of  uniaxial  strain.  In  order  to  avoid  am  energy¬ 
generating  loop  in  the  pres sure- volumetric  strain  relation,  the  constrained 
and  shear  mcuuli  must  be  restricted  by  the  following  inequality 


(3.8) 


To  avoid  negative  values  of  v  ,  the  value  of  G  is  further  restricted 
by  the  inequality 


2G  <  M* 


(3.9) 
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3.2  LINEAR  ELASTIC  G  MODEL  FOR  MATERIALS  WITH  A  NONLINEAR 

HYSTERETIC  M 

This  model  is  the  nonlinear  version  of  the  bilinear  hysteretic  model 
presented  in  Section  3.1.  Figures  3.2a  and  b  show  the  vertical  stress- 
strain  and  the  principal  stress  difference- strain  difference  relations, 
respectively,  for  this  model  under  conditions  of  uniaxial  strain.  The  ma¬ 
terial  first  behaves  as  a  linear  elastic  G  model  with  a  nonlinear  hys¬ 
teretic  constrained  modulus  (Section  2.5)  until  the  plastic  yield  stress 
is  reached.  This  corresponds  to  Point  1  in  Figure  3.2.  Upon  yielding, 
the  slope  of  the  principal  stress  difference- strain  difference  curve 
breaks  and,  becomes  equal  to  zero.  The  slope  of  the  vertical  stress- strain 
curve,  however,  increases  with  continued  application  of  vertical  stress. 

The  magnitude  of  the  stress  difference  at  Point  1  is  V3  k  ,  where  k  is 
the  von  Mises  yield  limit.  The  magnitude  of  the  vertical  stress  at  Point  1 
is  the  value  of  a  corresponding  to  e  =  ^3;  ♦  During  unloading  from 

a  stress  state  corresponding  to  Point-  2,  the  material  behaves  as  a  linear 
elastic  solid  in  shear  and  exhibits  nonlinear  hysteretic  behavior  in  the 
hydrostatic  part  of  the  deformation,  until  the  lower  yield  surface  corre¬ 
sponding  to  Point  3  is  reached.  The  slope  of  the  stress  difference- strain 
difference  curve  breaks  again  at  Point  3  and  becomes  equal  to  zero.  The 
slope  of  the  vertical  stress- strain  curve  continuously  decreases  with  fur¬ 
ther  unloading  of  the  vertical  stress  beyond  Point  3. 

This  description  of  the  vertical  stress- strain  relation  is  more  real¬ 
istic  than  the  bilinear  representation  presented  in  Section  3.1,  Since 
the  constrained  modulus  in  this  case  depends  on  the  stress  level,  the  con¬ 
struction  of  the  model  must  be  carried  out  incrementally  utilizing  Equa¬ 
tions  2.3  and  2.4  and  imposing  the  condition  that  (a  -  a  )  =  0  on  the 

z  r 

yield  surface. 

The  constructed  pressure-volumetric  strain  and  stress  difference- 
pressure  relations  for  the  model  are  plotted  schematically  in  Figures  3.2c 
and  d,  respectively.  In  order  for  the  pressure- volume trie  strain  curve 
(Figure  3.2c)  to  avoid  exhibiting  an  energy- generating  loop  and  to  unload 
along  Path  2-B,  the  following  inequality  must  hold  in  any  stress  increment 


(3-10) 


(MU  •  >  3  0 

As  was  pointed  out  previously,  the  above  inequality  insures  that  the  un¬ 
loading  slope  of  the  pres sure- volumetric  curve  is  greater  than  the  corre¬ 
sponding  loading  slope  in  any  stress  increment.  Ihe  shear  modulus  is 
further  restricted  to  positive  values  of  v  by  the  inequality 


2G  <  M. 
SL 


(3-11) 


in  any  stress  increment.  When  coupled  with  Equation  1.4,  the  above  in¬ 
equality  also  restricts  the  slopes  along  the  uniaxial  strain-stress  paths 
plotted  in  principal  stress  difference  versus  pressure  space  to  positive 
values  less  than  3>  i.e. 


az  "  °r 


<  3 


'  uniaxial 
strain 


(3.12) 


The  stress  path  in  Figure  3*2d  has  a  curvature  in  the  elastic  range 
(Path  0-1)  that  as  opposite  to  that  in  the  corresponding  portion  of  the 
vertical  stress-strain  curve  (Figure  3.2a).  Thus,  while  the  vertical 
stress- strain  curve  stiffens  as  pressure  increases,  the  stress  path 
softens  until  the  material  yields .  The  unloading  stress  path  between 
Points  2  and  3  for  this  model  appr  aches  a  straight  line  at  high  pres¬ 
sures  where  the  unloading  constrained  modulus  is  essentially  linear  and 
becomes  convex  to  the  pressure  axis  at  low  pressure  where  the  hook  or 
tail  in  the  vertical  stress- strain  curve  starts  developing.  Thus,  if  the 
constrained  modulus  is  continuously  hardening,  the  slope  of  the  calcu- 
daned  stress  path  is  continuously  reduced;  if  the  constrained  modulus 
softens,  the  calculated  stress  path  becomes  progressively  steeper.  This 
latter  calculated  curvature  when  unloading  reaches  low  pressure  is  in 
direct,  contrast  with  available  experimental  uniaxial  strain-stress  path 
data. 
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e.  PRESSURE-VOLUMETRIC  STRAIN  RELATION 


(<7«  "ar> 


b.  PRINCIPAL  STRESS  DIFFERENCE- 
STRAIN  DIFFERENCE  RELATION 


<o,-cr(> 


c.  PRINCIPAL  STRESS  DIFFERENCE- 
PRESSURE  RELATION  (STRESS 
PATH) 


Figure  3.1  Behavior  of  linear  elastic  shear  modulus  model  $jt 
elastic-plastic  materials  with  a  bilinear  hysteretic  constrained 
modulus  under  conditions  of  uniaxial  strain. 
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a.  VERTICAL  STRESS-STRAIN  RELATION 


P 


e  PRESSURE-VOLUMETRIC 
STRAIN  RELATION 


b.  PRINCIPAL  STRESS  DIFFERENCE- 
STRAIN  DIFFERENCE  RELATION 


(a,  -  a,) 


d.  PRINCIPAL  STRESS  DIFFERENCE- 
PRESSURE  RELATION  (STRESS 
PATH) 


Figure  3.2  Behavior  of  linear  elastic  shear  modulus  model  for 
elastic-plastic  materials  with  a  nonlinear  hysteretic  constrained 
modulus  under  conditions  of  uniaxial  strain. 
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CHAPTER  4 


CONCLUSIONS  AND  RECOMMENDATIONS 

4.1  CONCLUSIONS 

The  results  from  this  simple  examination  of  several  highly  idealized 
material  descriptions  indicate  that  many  of  the  behavior  characteristics 
of  real  earth  materials  often  recorded  during  laboratory  tests  can  be  sim¬ 
ulated,  at  least  qualitatively,  by  elastic-plastic  constitutive  models  of 
the  constant  G-type  by  observing  certain  restrictions  on  the  material 
property  parameters  G  and  M  .  Other  characteristic  phenomena,  however, 
cannot  be  mirrored  with  this  type  model  under  any  circumstances. 

For  example,  even  if  is  always  greater  than  M^  at  any  stress 
level,  only  by  utilizing  the  linear  hysteretic  G  model  and  observing 
the  following  inequalities, 

\  U  /  \  G  M 

(\-M*H(Gu-GJ  ^  G^ 

can  one  insure  that  pressure-volumetric  strain  curves  will  be  constructed 
that  exhibit  energy- absorbing  hysteretic  behavior  rather  than  energy¬ 
generating  loops  when  unloaded  from  any  stress  level .  These  same  inequal¬ 
ities  must  also  be  adhered  to  if  unloading  paths  of  principal  stress  dif¬ 
ference  versus  pressure  for  a  state  of  uniaxial  strain  are  to  fall  below 
the  loading  path  rather  than  above  it.  In  order  to  be  consistent  with 
the  experimental  observation  that  axial  and  radial  stresses  simultane¬ 
ously  increase  or  decrease  during  loading  or  unloading,  respectively,  and 
also  restrict  Poisson's  ratio  values  to  be  greater  than  zero  at  any  stress 
level,  it  is  necessary  to  observe  the  inequalities 

2G .  <  M.  and  2G  <  M 
l  &  u  u 

The  linear  hysteretic  G  model  will  obviously  incorporate  the  well  known 
ability  of  soils  to  dissipate  energy  in  shear,  even  at  subyield  stress 
levels . 

Uniaxial  strain  tests  on  undisturbed  natural  soils,  however,  usually 
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reveal  a  reversal  of  curvature  in  a  versus  e  plots  during  loading, 

2  2 

first  softening  and  then  hardening,  •while  the  corresponding  (a  -  a  )  ver- 
sus  p  plots  generally  indicate  monotonically  softening  behavior.  Any  con¬ 
stant  G  model,  however,  if  coupled  with  a  loading  constrained  modulus 
that  reverses  curvature,  will  of  necessity  have  a  curvature  reversal  during 
loading  in  its  (c  -  a  )  versus  p  curve,  unless  a  plastic  yield  criterion 

2  i* 

is  invoked  at  a  pressure  level  less  than  that  of  the  inflection  point.  On 
the  other  hand,  reversals  of  curvature  have  never  been  observed  in  plots  of 
unloading  constrained  moduli;  yet  unloading  stress  paths  frequently  exhibit 
directional  changes  in  curvature.  Again,  if  the  constrained  modulus  'un¬ 
loading  curvature  does  not  change  sign,  then  the  unloading  stress  path  cur¬ 
vature  calculated  with  any  constant  G  model  cannot.  In  fact,  when  unload¬ 
ing  reaches  low  pressure  or  approaches  the  lower-bound  or  negative  plastic 
yield  limit,  the  constant  G  calculated  stress  path  curvature  is  always  in 
direct  contrast  with  available  experimental  uniaxial  strain  test  data- 

While  it  is  obvious  that  elastic-plastic  constitutive  models  formu¬ 
lated  with  constant  values  of  G  can  match  some  of  the  dominant  features 
exhibited  by  soils  and  weak  rocks  during  uniaxial  strain  tests,  it  is 
equally  obvious  that  some  feature's  can  only  be  mirrored  in  a  limited  fash¬ 
ion  and  others  not  at  all.  It  is  also  quite  obvious  that  any  model  that  is 
linear  in  G  cannot  possibly  exhibit  the  nonlinear  stress- strain  behavior 
usually  observed  during  triaxial  compression  tests. 

4 . 2  RECOMMENDATIONS 

Thus,  it  is  recommended  that  research  efforts  to  develop  realistic 
mathematical  constitutive  models  of  earth  media  for  use  in  free- field 
ground  shock  calculation  codes  continue.  In  the  authors'  opinion,  such  ef¬ 
forts  should  be  governed  by  the  following  criteria: 

1.  The  models  should  be  capable  of  qualitatively  and  quantitatively 
matching  the  salient  nonlinear  and  hysteretic  response  characteristics  of 
earth  media,  not  only  us  determined  from  the  uniaxial  strain  test,  but  also 
as  determined  from  a  variety  of  other  laboratory  test  boundary  conditions. 
If  confidence  is  to  be  placed  in  a  ground  shock  code  which  must  calculate 
earth  material  response  to  an  unknown  number  of  varying  states  of  stress, 
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transient  boundary  loadings,  and  erratic  site  geometries,  the  material 
model  used  should  at  least  have  the  capability  to  reasonably  duplicate  the 
response  of  relatively  '  .Imogen  eous  specimens  of  simplified  geometries 
loaded  in  the  laboratory  along  known  stress  paths  under  carefully  controlled 
conditions . 

2.  The  numerical  values  of  the  coefficients  defining  the  models 
should  be  readily  derivable  from  laboratory  test  data.  It  is  most  impor¬ 
tant  that  the  coefficient  values  not  be  merely  a  set  of  numbers  generated 
through  a  trial  and  error  "black  box"  routine  to  fit  a  given  set  of  data, 
but  that  they  have  physical  significance  in  terms  of  compressibility,  shear 
strength,  etc.,  so  that  when  extrapolating  to  different  materials,  rational 
engineering  judgments  can  be  made  as  to  their  relative  magnitudes  based  on 
geologic  descriptions,  mechanical  properties,  and  other  conventional 
indexes . 

3.  Theoretical  restrictions,  such  as  those  imposed  by  uniqueness  and 
continuity  considerations,  should  be  satisfied  to  the  maximum  practicable 
extent.  It  should  be  kept  in  mind,  however,  that  earth  materials  are  al¬ 
most  invariably  nonhomogeneous ,  anisotropic,  stress-history  dependent, 
chemically  entangled  mixtures  of  gases,  liquids,  and  solid  particles. 
Therefore,  it  seems  unreasonable  to  expect  that  constitutive  models  for 
these  materials  can  be  developed  which  satisfy  all  the  mathematical  re¬ 
strictions  imposed  through  assumption  of  formal  continuous  media  theories 
while  accurately  mirroring  all  of  the  details  observed  during  various  lab¬ 
oratory  tests.  Obviously,  some  compromises  will  inevitably  have  to  be 
made;  parametric  studies  to  evaluate  the  relative  importance  of  ignoring 
certain  features  of  either  data  or  theory  are  needed. 

4.  Any  ground  shock  prediction  technique  or  code  should  be  thoroughly 
checked  against  the  available  analytical  solutions  and  a  variety  of  actual 
field  experiments.  Efforts  to  incorporate  new  constitutive  models  into 
generalized  codes  for  calculating  transient  ground  shock  effects  from  ex¬ 
plosive  events  should  also  include  attempts  to  obtain  analytical  or  closed- 
foim  solutions  to  special  problems  of  wave  propagation  through  materials 
containing  the  principal  features  of  the  new  models .  Presumably,  improved 
models  should  provide  improved  correlations  with  existing  field  data. 
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'ihe  above  criteria  are  admittedly  more  pragmatic  than  theoretical. 
However,  the  authors  strongly  believe  that  any  critical  analysis  of 
material  models  should  include  examination  in  the  light  of  all  these 
criteria  in  order  to  maintain  a  balance  in  perspective. 


APPENDIX  A 


CONSTITUTIVE  EQUATIONS  AND  UNIAXIAL  STRAIN  BEHAVIOR 
FOR  A  PRANDTL- REUSS  MATERIAL 


A.l  PURPOSE 

When  constructing  constitutive  models  to  represent  the  rather  compli¬ 
cated  nonlinear  hysteretic  behavior  of  real  earth  materials  using  incre¬ 
mental  stress-strain  relations  based  on  classical  elastic-plastic  theory, 
it  is  helpful  to  review  the  equations  governing  the  behavior  of  an  ideal 
elastic-plastic  material  and  examine  its  response  under  typical  test 
boundary  conditions. 

The  purpose  of  this  appendix  is  to  perform  such  a  function  for  per¬ 
haps  the  most  widely  known  ideal  elastic-plastic  material,  i.e.  the 
Prandtl-Reuss  material,  and  the  most  useful  test  conditions  for  blast- 
oriented  problems,  i.e.  loading  and  unloading  states  of  uniaxial  strain. 

A. 2  CONSTITUTIVE  EQUATIONS 

A. 2.1  Yield  Criterion.  The  yield  criterion  associated  with  the 
Prandtl-Reuss  material  is  the  well-known  von  Mises  criterion  given  by 


where  k  is  a  material  constant  related  to  the  cohesive  strength  of  the 
material.  Equation  A.l  describes  a  right- circular  cylinder  in  the  prin¬ 
cipal  stress  space  with  its  central  axis  the  line  of  hydrostatic  stress 
as  shown  in  Figure  A.l. 

According  to  the  basic  postulate  of  the  plasticity  theory,  the  mater 
rial  would  flow  plastically,  undergoing  plastic  as  well  as  elastic  strains 
when  Equation  A.l  is  satisfied.  The  total  strain  is  then  the  sum  of  the 
elastic  (recoverable)  strain  and  plastic  or  permanent  strain.  When  the 
stresses  are  less  than  those  satisfying  Equation  A.l,  the  material  will 
undergo  elastic  strains  only;  i.e.,  the  material  behaves  as  a  linear 
elastic  solid.  In  incremental  notation,  these  conditions  can  be  written 
as 
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•  •  e 

e.  .  =  er . 

ij  13 


when 


1.1  4  id 


\FI<k 

when  >Fk =  k 


(A.2) 


:e 


whe re  e.  .  =  components  of  elastic  strain  increment  tensor 

3.J 


e?. 

13 


components  of  plastic  strain  Increment  tensor 
According  to  Hookes'  law,  the  elastic  strain  increment  tensor  is  given  by 


ee  -  3;  *  v  *  v  *  R 

ij  E  °ij  “  E  °kk  ij 


(A. 3) 


where  E  =  Young's  modulus  of  elasticity. 

The  plastic  strain  increment  tensor  can  be  derived  from  the  yield 
criterion  (Equation  A.l)  by  utilizing  the  theory  of  plastic  potential: 


P  =  x-M. 
id  <*.. 


(A.4) 


where  \  =  a  positive  scalar  factor  of  proportionality 
/  =  J  ’  -  k2 

A. 2. 2  Stress-Strain  Relationship.  In  view  of  Equation  A.4  and  the 
function  /  ,  the  plastic  strain  increment  tensor  becomes 


'k  - k2) 


=  X 


=  Xs 


id 


(A. 5) 


Since  s^  =  0  ,  it  follows  from  Equation  A. 5  that 

«*-  0 
kk 


(A. 6) 


indicating  that  no  plastic  volume  change  can  occur  in  the  plastic  range 
for  Prandtl-Reuss  material.  Hence,  all  volume  change  is  elastic,  i.e. 


'kk  “  'kk 


(A. 7) 
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and 


p  r:  K 6 


:e 

kk 


(A.8) 


Accordingly,  the  plastic  deviator  strain  increment  tensor  et\  is  the 

same  as  the  plastic  strain  increment  tensor  6? .  .  Ihe  total  deviator 

strain  increment  tensor  e.  .  then  becomes 

10 


•  *e  .  *p 
e.  .  -  e.  .  +  £ .  . 
10  10  10 


(A.9) 


•  0  •  ■ 

where  e .  .  =  s.  ./2G  =  components  of  elastic  deviator  strain  increment  ten- 
ij  10 

sor.  Employing  Equation  A. 5  in  Equation  A.9  results  in  the  following  re¬ 
lation  for  the  total  deviator  strain  increment  tensor 


e  =  fii  +  Xs 

6ij  2G  XSij 


(A.10) 


Multiplying  both  sides  of  Equation  A. 10  by  2Gs. .  yields 

10 


2Gs_.  ^e,.  4  =  s_.  ..s..  4  +  X^Gs4,s 


j  •  •  w  •  »  u  >  •  o  •  • 

1J  10  10  10 


iO  i0 


(A. 11) 


or 


2GW  =  +  4G\Jg 


(A. 12) 


where  W  -  s^e^j  =  the  rate  of  work  or  increment  of  internal  energy  due 

to  distortion  and  Ji  =  s.  .s,  .  —  Tt  -  1-2 

2  io  id 

ing,  Jg  =  0  and  Equation  A. 12  yields  the  following  relation  for  the  pro¬ 
portionality  factor  X 


Since  Jg  =  k  =  constant  during  yield- 


X  =  A  (A.13) 
2k 

It  is  to  be  noted  that  W  is  only  due  to  plastic  deformation.  In 
view  of  Equation  A.13,  the  total  deviator  strain  increment  tensor  (Equa¬ 
tion  A.10)  takes  the  following  form 
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(a.i4) 


e.  . 
10 


Equation  A.l4  can  be  inverted  to  solve  for  the  deviator  stress  increment 

tensor  s.  .  ,  i.e. 

10 


s.  . 
10 


=  2G  e.  . 

V  10 


(A.15) 


Equation  A.15  (or  Equation  A.l4)  is  the  Prandtl-Reuss  stress-strain  rela- 

2 

tion;  it  applies  in  the  plastic  range  when  Jg  =  k  and  W  >  0  .  If 
W  <  0  ,  unloading  is  taking  place,  i.e.  the  behavior  is  only  elastic  and 
Equation  A. 3  applies.  If  W  =  0  the  loading  is  said  to  be  neutral; 

The  above  results  can  be  summarized  in  the  following  form 


e-.,  =  -M-  when  J J*  <k  (elastic  (A.l6a) 

*•  loading) 

hz  +  f  sij)  when'  =  k  »  W  >  0  (A.i6b) 

*  ^3  i  •  +  '- 

e.  .  =  -rpf  when  W  <  0  (elastic  (A.l6c) 

J  ^  unloading) 


In  order  to  obtain  the  components  of  the  stress  increment  tensor 
Equation  A.l6  must  be  combined  with  Equation  A. 8. 


A. 3  BEHAVIOR  UNDER  CONDITIONS  OF  UNIAXIAL  STRAIN 


For  uniaxial  strain  conditions  in  a  cylindrical  coordinate  system 
r,0.z  ,  the  following  boundary  conditions  prevail  (Figure  A. 2) 


e  =  e  =  o  > 

r  9 


e 

r 


(A.17) 


e 

z  3  z  J 

In  the  elastic  range,  the  behavior  of  the  material  is  governed  by  Equa¬ 
tion  A. 3,  which  can  also  be  written  in  the  following  form 
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(A.18) 


€.  , 
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s.  . 
2G 


JL 

3k 


S.  . 

3-0 


•  < 

Since  £  =  e  =  0  ,  Equation  A.l8  yields  the  following  relations  between 
r  0  * 

components  of  deviator  stress  increment  and  p 


s 

r 


In  view  of  the  fact  that 


Equation  A. 19  gives 


(A.19) 


(A. 20) 


•  •  ♦ 

Utilizing  the  relation  s  =  a  -  P  in  Equation  A. 20  results  in  the  fol- 

z  z 

lowing  relation  for  the  vertical  stress  increment  a 

% 

s  ■  (*£*>  <*•*> 

Since  for  uniaxial  strain  conditions  ,  and  in  view  of  Equations 

A. 7  and  A. 8,  the  vertical  stress-strain  increment  relation  in  the  elastic 
range  becomes 


(A. 22) 


From  Equation  A.19,  the  principal  stress  difference,  in  incremental  form, 
is 


(•  •  \  or*  • 

°z  ■  °r/  T  P  (A. 23a) 

indicating  that  the  slope  of  the  uniaxial  strain-stress  path  in  the 
elastic  range  is  equal  to  2Ci/k  .  From  the  relations  between  constants 
in  elastic  theory,  it  can  readily  be  shown  that 

f  ■  (*•**) 
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which ,  when  coupled  with  Equation  A  .23a,  gives  the  following  relation  for 
the  condition  of  uniaxial  strain 


(A. 24) 


In  view  of  Equations  A. 8,  A. 17,  and  A. 23,  the  principal  stress  difference- 
strain  difference  relation  for  uniaxial,  strain  becomes 

(az  -  or)  =  2G(ez  -  =  20^  (A. 25) 

For  virgin  loading  in  the  elastic  range,  Equations  A. 17  through  A. 25 
can  be  used  without  the  dot  notation. 

In  the  uniaxial  strain  test 


n/jJ  (®z  -  or) 


V3 


Thus  the  material  wil],  yield  when 


(A. 26) 


*)-k 


(A. 27 ) 


In  view  of  Equations  2.8,  A. 24,  and  A. 27,  the  value  of  vertical  stress 
cv  at.  yield  becomes 


(A. 28) 


Thus,  when  a  reaches  the  value  given  by  Equation  A. 28,  the  material 
z 

yields  and  the  continued  application  of  vertical  stress  causes  the  mate¬ 
rial  to  move  along  the  yield  surface  undergoing  plastic  as  well  as  elas¬ 
tic  strains.  In  the  plastic  range  Equations  A.l8  through  A. 25  no  longer 
apply  and  recourse  to  Equations  A. 8  and  A.15  is  necessary. 

« 

According  to  Equation  A.15,  the  deviator  stress  increment  s  in 

z 

the  plastic  range  is  given  by 
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(A. 29) 


*  GW 

S  ~  2Ge - -  s 

z  z  2  z 

K 

The  rate  of  work  W  for  conditions  of  uniaxial  strain  reduces  to 


W  =  s  e  +  2s  e 
z  z  r  r 


(A. 30) 


or,  utilizing  the  fact  that  s^  =  c^  =  0  , 


W  =  |  s  e 
2  z  z 


(A. 31) 


Substituting  Equation  A. 31  into  Equation  A. 29  and  utilizing  the  relation 

•  2  •  * 
e  =  —  e  results  in  the  following  expression  for  s 
Z  j  z  z 


s  ~ 


4g 
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2* 
s  e 

z  z 


(A.32) 


In  the  plastic  range 


k2  = 


J2 


_  1  / 

“3  (as“ 
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(A. 33) 


and  Equation  A, 32  reduces  to 


s 

z 


4 

3 


z 


=  0 


Since 


s  =  0  also,  and 
r  ' 


a 

z 


sz  +  P  =  P 


(A. 34) 


(A. 35a) 
(A.  35b) 


which  indicates  that  the  material  behaves  as  though  it  were  a  fluid  once 
it  has  reached  its  limiting  shear  resistance.  Substituting  Equation  A. 8 
into  Equation  A. 35a  and  considering  the  fact  that  e  =  e  ^  the  verti- 
cal  stress-strain  increment  relation  in  the  plastic  range  becomes 
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(A. 36) 


Thus  ,  the  loading  slope  of  the  a  versus  e  curve  breaks  or  softens 

2  z 

when  yielding  occurs  and  becomes  equal  to  the  bulk  modulus.  Accordingly, 
the  loading  slopes  of  the  principal  stress  difference-pressure  curve  and 
the  principal  stress  difference- strain  difference  curve  become  zero.  The 
slope  of  the  pressure-volumetric  strain  curve,  however,  is  not  affected 
and  ranains  constant. 

Once  the  material  unloads,  it  behaves  as  a  linear  elastic  solid 
again,  satisfying  Equations  A.l8  through  A. 25.  If  unloading  is  contin¬ 
ued  until  the  lower  yield  surface  corresponding  to 


-fc(az 
V3  2 


or)  =  -k 


(A. 37) 


is  reached,  the  material  flows  plastically  again  and  Equations  A. 8  and 
A. 15  govern  the  behavior  of  the  material. 

The  foregoing  analyses  are  depicted  schematically  in  Figure  A. 3. 

From  Figure  A. 3  it  can  readily  be  seen  that  a  Prandtl-Reuss  material  ex¬ 
hibits  no  hysteretic  effects  during  a  hydrostatic  test.  Furthermore,  for 
a  Prandtl-Reuss  material,  the  vertical  stress-strain  curve  from  a  uniax¬ 
ial  strain  test  would  break  or  soften  when  yielding  occurs  and  remain 
concave  to  the  strain  axis  with  continued  application  of  vertical  stress. 
These  observations  are  generally  contrary  to  existing  experimental  data 
on  the  stress- strain  behavior  of  earth  materials  even  from  a  qualitative 
point  of  view.  This  simply  illustrates  a  fact,  of  no  great  surprise,  i.e. 
that  the  ideal  linear  elastic-plastic  (Prandtl-Reuss)  model  does  not  ade¬ 
quately  describe  the  stress- strain  behavior  of  real,  nonideal  soils — in 
particular  the  behavior  under  conditions  of  uniaxial  strain. 
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Figure  A. 2  Uniaxial  strain  configuration  in 
cylindrical  coordinate. 


o.  VERTICAL  STRESS-STRAIN  RELATION 


c.  PRESSURE-VOLUMETRIC 
STRAIN  RELATION 


b.„  PRINCIPAL  STRESS  DIFFERENC 
STRAIN  DIFFERENCE  RELATION 


d.  PRINCIPAL  STRESS  DIFFERENCE- 
PRESSURE  RELATION  (STRESS 
PATH) 


Figure  A.  3  Behavior  of  Prandtl-Reuss  material  under 
conditions  of  uniaxial  strain. 
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